The Importance of Intergalactic Structure to Gravitationally Lensed Quasars 
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ABSTRACT 

■ Image flux ratio anomalies have been attributed to substructures within the gravitational 
I lens and to small mass halos (M ^ 10^° Mq) in intergalactic space. In this paper, analytic 

calculations are presented that help in the understanding of how intergalactic halos affect 
magnification ratios. It is found that intergalactic halos can produce anomalies at a similar 
^ [ level to those that are observed. Intergalactic halos with masses < 10^° M0 are expected to 

' cause relative deflections between images of order 10 milliarcseconds, which are then magnified 

, by the primary lens. They will also cause fluctuations in the surface density on the several 

CSJ I percent level. The importance of intergalactic halos depends strongly on the radial profile of the 

halos and the primordial power spectrum at small scales. Strongly lensed quasars provide an 
^ , opportunity to probe these properties. A strong dependence on the QSO rcdshift is predicted 

00 ' and can be used to distinguish between intergalactic structure and substructure as the cause of 

' magnification anomalies. This analytic approach also explains why some previous semi-analytic 

, estimates disagreed with numerical calculations. 

■ 1. Introduction 

i-C ' Gravitationally lensed quasars have recently been used to identify substructure in the dark 



matter halos of lens galaxies (Mao & Schneider 1998; Metcalf & Madau 2001; Metcalf & Zhao 



O ■ 2002; Metcalf 2002; Chiba 2002; Dalai & Kochanek 2002; Bradac et al. 2002; Keeton 2003; 

' Metcalf et al. 2004). The Cold Dark Matter (CDM) model predicts a certain amount of small 

. mass subhalos ( ;$ 10*" Mq) that do not have visible dwarf galaxy counterparts. These subhalos 

within the gravitational lens could be causing the observed anomalies in the image flux ratios, 
but it is also possible that they are being caused by clumps of matter somewhere along the line 
of sight, but not associated with the lens galaxies. Mao et al. (2004) have argued that ACDM 
' does not predict enough substructure within galactic halos to account for the anomalies. Amara 

et al. (2004) have studied the lensing properties of lenses taken out of CDM simulations and 
found that they do not contain enough substructure, but this conclusion may be dependent on 
resolution. Future cosmological simulation with higher resolution might be able to verify this 
for certain. 

It was found in Metcalf (2004), using numerical simulations, that intergalactic halos could 
account for most of the anomalies seen in the magnification ratios. This conclusion has been 
somewhat controversial largely because it disagrees with the conclusions of Chen, Kravtsov, & 
Keeton (2003) who found that the contribution from intergalactic halos would be ^ 10% of 
that from substructure within the primary lens halo, not enough to account for the observed 
anomalies by themselves. Chen et al. (2003) relied on an analytic calculation where the 
deflection caused by the primary lens is approximated by a first order expansion and then the 
cross section as a function of magniflcation is calculated for a single perturbing substructure. In 



1 Hubble Fellow 



- 2 - 



this paper, it is demonstrated why this type of calculation is inadequate and how the results of 
the numerical simulations (Metcalf 2004) can be understood in terms of a simple analytic model. 

This paper also makes clear what properties of small scale structure most affect the 
magnification ratios of strongly lensed quasars. Previously, the perturbing halos were simply 
modeled as truncated Singular Isothermal Spheres (SISs), except in Metcalf (2004), while here 
more realistic profiles are used and the sensitivity of lensing to their structure is investigated. 
The analytic model developed here demonstrates how flux ratio anomalies and Spectroscopic 
Gravitational Lensing (SGL) (see Metcalf et al. (2004)) can be used to study the small scale, or 
small mass, end of structure formation. 

For the calculations in this paper the standard ACDM cosmological model is used with 
Clm = 0.3, Oa = 0.7, Ho = 70 km s~^ Mpc~^ and a scale free initial power spectrum. In the 
next section the density profiles and number densities of small mass dark matter halos are 
discussed. In section 3 the importance of intergalactic structure to the magnifications of QSO 
images is calculated both using a halo model, section 3.1, and directly from the power spectrum 
of density fluctuations, section 3.2. The conclusions are discussed in section 4. 

2. Halo Profiles and Concentrations 

Prom cosmological N-body simulations it has been found that the density profile of a dark 

matter halo tends to take a universal form given by (at least approximately) Navarro, Frenk 
& White (NEW) - p{r) = ^(i+r'/r' )^ where pc is the critical density of the universe and is 
the scale size (Navarro, Frenk, & White 1997). The concentration of the halo is defined as 
c = r^oo/rs where r2oo is the radius within which the average density is 200 times the mean 
density of the universe, a proxy for the virial radius of the halo. The outer boundary of the halo 
is generally taken to be r2oo as will be the case here. Once c and the mass, M, of a halo are 
fixed, 5c and Vs can be calculated. If a relationship between c and M is specified, the profile is a 
one parameter family. 

The concentration increases with decreasing halo mass in the CDM model of structure 
formation. This is a result of smaller halos generally forming earlier, when the universe was 
denser. Navarro et al. (1997) propose a relation between concentration and mass by assuming 
that when a halo forms it contracts by a fixed factor from the mean density of the universe at 
that time. There is some ambiguity in defining the collapse time and not all halos of a given 
mass form at the same time so this proscription is at best approximate. More recently, fits to 
N-body simulations give approximately a power-law dependence of the concentration on subhalo 
mass. 



Values for /3 range firom -0.1 to -0.2 and ~ 14 (Bullock et al. 2001; Colin et al. 2004). In 

Metcalf (2004) a redshift independent c-M relation was used as in Navarro et al. (1997). N-body 
simulations do show that the concentration of pure dark matter halos of a fixed mass should 
have a redshift dependence given above. This decreases their strength as lenses somewhat, as 
will be seen in section 3. However, there are a number of factors the complicate the c-M relation 
and would be expected to increase the effective concentration of a halo. 

In addition to dark matter, halos contain baryons that can significantly change their profiles. 
For example, the contraction of the baryons because of cooling and angular momentum loss 
will increase their density at the center of the halo and drag in additional dark matter, making 
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the profile steeper and more concentrated. Lensing results indicate that the inner regions of 
large galaxies have a profile that approximates a singular isothermal sphere (SIS), p{r) oc r^^ 
(Treu & Koopmans 2004). These effects are expected to be smaller for small mass halos where 
the baryons do not cool efficiently and can be expelled by supernova feedback (Dekel & Silk 
1986). Observed dwarf galaxies are known to have higher mass to light ratios than L» galaxies 
(Mateo 1998; Kleyna et al. 2002). This also indicates that the profiles of small mass halos are 
more directly related to early structure formation driven by dark matter and less confused by 
feedback from star formation. With this in mind, it seems reasonable to approximate the profile 
of small mass halos by the NFW profile with a c-M relation approximating (1). but, given the 
uncertainties and for the purposes of comparison with Metcalf (2004), a redshift independent 
c-M relation will also be used in calculations. 

The c-M relation in N-body simulations also has a large scatter about the median (Bullock 
et al. (2001) find a 1 cr scatter of A(logc) ~ 0.2). This is likely to increase the importance of 
these halos for lensing since their effects generally increase with concentration more rapidly 
than linear. This is another reason why expression (1) is a conservative approximation for our 
purposes. 

The mass function of intergalactic halos is calculated using the Sheth-Tormen (Sheth & 
Tormen 2002) modification of the Press-Schechter (Press & Schechter 1974) method. This is 
known to agree with N-body simulations down to M ~ 10^ (Colin et al. 2004). This 

prescription accounts only for isolated halos. All halos that have been resolved in N-body 
simulations contain substructure within them. This hierarchy of subhalos could conceivably 
increase the importance of intergalactic halos to lensing. This will not be taken into account in 
this paper. 



3. Lensing by Intergalactic Halos 

The situation of interest is one where a large primary lens that is creating multiple images 
of the QSO and additional small perturbing halos along the line of sight. First we will use the 
multiple lens planes approximation (see figure 1). Here the deflections caused by each object 
are treated as if they take place suddenly in the plane of the lens and the light follows an 
unperturbed geodesic between the planes. This is known to be a very good approximation. 
Given the angular position of a point on the source, /3, we must calculate the image points, 0, 
that correspond to it. If there are N lens planes these angular positions are related by 

$,^, = ^e-j2^^d,(p,$) = $N+m (2) 

where Di is the angular size distance to the ith lens, Di,j is the distance between the ith and 
the jih lens planes and Dg — D^+i is the distance to the source. The position in physical units 
on the ith lens plane is ri = Ds<i>i and the deflection angle caused by that lens is di{fi). In 
general, the sub- or superscript I will refer to the primary lens. The distance to the primary lens 
is Di and Di,m+i = Dig. Equation (2) is only valid for a geometrically flat cosmology because it 
assumes that £),;,i+i +-Di+i,i+2 = -0^,^+2- Wc assume that this is the case in this paper. 

The magnification matrix is conveniently decomposed into the convergence, k, and two 
components of shear, 71 and 72, 

AUff\- _ f 72 \ , . 
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true light paths / = / 




Fig. 1. — A schematic diagram of the type of lensing system being considered. There is one primary lens that is responsible for 
the multiple images of the source. In addition, there are many secondary lenses (most not shown). The two images represented 
here are separated by an angle 9. There may also be two other images that are not shown. The unperturbed light paths are 
deflected only by the primary lens and with an appropriate model for the primary lens will meet on the source plane. If the 
deflections from secondary lens planes are taken into account without changing the primary lens model the light will follow 
the perturbed light paths (dotted curves) and will not necessarily meet at a single source. On the source plane they will be 
separated by a distance AaDg. To get the true light paths the primary lens model needs to be adjusted so that they again 
meet at a single point on the source plane. The angle Ao is thus a measure of how wrong the lens model is when intergalactic 
lensing is not taken into account. This diagram is not to scale in any respect. 

When there are multiple lens planes it is possible to get a rotation , but it will be small and 
of a higher order in the perturbations than is considered in this paper. In addition, rotation 
will not affect the magnification ratios. If there is only one lens plane, the convergence can be 
directly related to the local surface density, K{fi) — S(n)/Scrit- The critical surface density is 
5^crit = {4:TrGDiDis/c^Ds)^^ ■ The magnification of a point-like image is /i = (dct^)^^. 

3.1. simplified model 

The influence of intergalactic halos on the magnification ratios and the cusp caustic relation 
was studied in Metcalf (2004) by directly solving equations (2) in numerical simulations. 
Numerical simulations are necessary for making accurate quantitative predictions, but a 
simplified model can help greatly with our understanding of which parameters and effects are 
important. 
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The essential difficulty with solving the lens equation (2) is that the positions on each lens 
plane are dependent on the positions on the other lens planes. The strategy taken here is to 
calculate quantities along unperturbed lines of sight, or the path the light would take if only the 
primary lens were present. For an unperturbed path, the position on the ith plane is 

; w 

This approximation is valid as long as the deflections by intergalactic halos are small so that 
quantities that are second order in their deflections and derivatives of their deflections can be 
ignored. 

The surface density of the NFW profile can be calculated directly, but the result is 
cumbersome. A simplified surface density can be used that exhibits the essential properties of 
the NFW profile and allows for some adjustments 

S(r,2) = -— (5) 

where n ~ 2 is most like the NFW profile although n = 3 will also be used to gain some insight 
into how the profile affects the lensing properties. For n = 2 the mass within the radius r is 



M{x) = 2-KrlT.o 



In (1 + x) - 



l + x 



T 

(6) 



The total mass is -M200 = M{c) where the concentration, c(M2oo), is given by (1). For n = 3 a 
different expression for M{x) is derived, but M200 and c(Af2oo) arc defined in the same way. 

Small halos affect the magnification of an image in two ways. One is from the halos changing 
the image positions and the primary lens' magnification being different at the new locations. 
Another way is for the substructure to contribute directly to the magnification of the image 
by changing A^^iO) without changing the image position enough to change the host lens's 
magnification. One of the goals of this paper is to determine which of these effects is more 
important. For this purpose the deflection, the convergence and the shear along the unperturbed 
path will be considered separately in the next two subsections. 



3.1.1. deflection 

The deflection angle caused by a axisymmetric halo is 

a{r) = AGM{r/rs)^. (7) 

From the lens equation (2) it is seen that the important quantity is d(r) = ^^^j^Q.{r) which is 
the lowest order displacement of a light ray on the source plane (see figure 1). 

Since \ct{r)\ falls only as r~^, for large r it might be expected that the deflections of many 
halos will add up to a total that is substantially larger than what an individual halo is likely to 
cause by itself. In fact, {a{6)'^) diverges logarithmically in the plane lens approximation because 
the number of halos rises as r^. This is not a problem here however because only differences in 
the deflection angle are important. A constant deflection over the entire lens is equivalent to a 
shift in the position of the source. Such a uniform shift would be taken into account when fitting 
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Fig. 2. — The root-mcan-squared perturbation to the angular separation, Ad((?), as a function of images separation, 9. The 
range of halo masses included is 10'* Mq < M < IQi^^ M© for the upper set of curves, 10* M© < M < 10^° Mq for the middle 
set and 10'* Mq < M < 10® Mq for the lowest set. In each set there are curves with varying values for the power spectrum 
normalization, as, and the dependence of concentration on halo mass, /3. 



a primary lens model so it eannot produce any anomaly in the magnification ratios. Howc;vcr, if 
the deflections are different for the different images then they will cause the image positions to 
shift independently. 

To asses the significance of these deflections we calculate the correlation function between 

the perturbations to the deflection caused by intcrgalactic halos as a function of the image 
separation. The positions of the halos will be treated as random and uncorrelated. In this case 
the correlation is 



Aa{ef = ( a{0') - a{0' + 0) 



_ 1 r" 



dz 



D{z,z,f 



^/Q{1 + z)^ + 
j dPr \d{r) - d{r + y{z,e))\^ 



Jo 



dN 



(8) 



(9) 



where y(z, 9) is the physical separation of two unperturbed lines of sight given by (4). These 
integrals are calculated numerically. The range of the r integral in (9) is increased until the 
answer stops changing. 

Figure 2 shows Ad{0) for a source at = 3 and lens at zi = 0.5. The most important 
parameter for Aq!(^) is the upper halo mass cutoff, Mmax- The concentration and ag are 
relatively unimportant. Using a redshift independent c-M relation gives slightly larger values 
for Aq!(^), but it is much smaller than the dependence on Mmax- Figure 3 what redshifts are 
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Fig. 3. — The contribution to Aa(0)^ from intcrgalactic hales as a function of redshift. The curves that end at z = 3 are 
for a source at that redshift and a primary lens at zi = 0.5. For the others Zs = 2 and 2; = 0.3. The halo mass range is 
10* Mq < M < lO'^'^ Mq in all cases except for the short dashed curves where 10* Mq < M < 10^ Mq. The image separation 
is set to 1 arcsecond. 



contributing to Aa{6). Four a QSO at ^ = 3 the majority of the deflection is coming from halos 
above z = 1. 

The deflections shown in figure 2 can have a significant effect on the image magnifications 
even if the convergence and shear caused by the substructures are small. This is because, for a 
single image, the perturbation to the deflection acts like a local change in the source position. 
The change in the local source position can be magnified by the host lens to resulting in a 
substantial change in the image position. Even if the intcrgalactic halos did not affect the local 
magnification as a function of image position they would affect where the images form and thus 
change the magnification caused by the primary lens. This is an effect that is completely absent 
from a cross section type calculation where the host lens' deflection is expanded to flrst order 
(constant magnification). Local changes in the source position on the order of 15 milliarcseconds 
do cause changes in the magnification ratios even in Einstein cross type lenses. However, as will 
be shown in section 3.1.2, the intcrgalactic halos will modify the magnification directly as well 
and probably this second effect is more important. 

The clustering of the halos has not been taken into account in this calculation. Clustering 
into filaments and other structures will only increase the value of Aa{6)'^. Because the deflection 
caused by a halo has such a long range the increase in Aa{9)^ could be signiflcant, but it will 
not be investigated here. 
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Fig. 4. — The standard deviation in the dimensionless surface density in small scale structure along random lines of sight as 
a function of source redshift. For the right panel the power spectrum normalization is erg = 0.9 as measured by (Spergel et al. 
2003) and on the left a lower normalization is used as = 0.7. The halo mass-concentration relation is specified by the /3 values 
for each type of curve. For each set of like curves the upper one is for a halo mass range of 10^ Mq < M < 10^^ Mq, the 
middle one is for 103 M© < M < 10i° Mg and the lower one is for 10^ M© < M < 10^ Mq. 



3.1.2. convergence & shear 

By taking an unperturbed light path in the lens equation (2) and the then calculating its 
derivatives it is easily shown that, to lowest order, the total convergence, k, is the sum of all the 
surface densities of all the halos along the line of sight weighted by the critical density. 

In the standard halo model of structure in the universe, all mass is in halos of some size. The 
standard expression for the angular size distance in a homogeneous universe takes into account 
the average surface density so, what is important for perturbations in the image magnifications 
is the variation in k,{9) about the mean. In addition, if the perturbation is uniform over the 
entire lens area, it will not affect any observables because of the mass sheet degeneracy. To 
cause anomalies in the magnification ratios k{6) must have significant fluctuations on a scale 
that is smaller than rotighly the Einstein ring radius of the primary lens 1 arcsecond). 

From the halo profile and the mass function the variance in n along an unperturbed line of 
sight, af. = {k^) — (k)^, can be calculated assuming uncorrelated halo positions in the same 
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Fig. 5. — The standard deviation in the dimensionless surface density in small scale structure. Everything is the same as 
in figure 4 except that the halo profile is modified so that it is cutoflF more sharply, n = 3 in equation (5). The variance is 
significantly larger. 

manner as the deflection angle was calculated in section 3.1.1, 



The mass integral is boimdcd on the top by the minimum mass of a halo that one would expect 
to have an observable galaxy inside it. Any observed galaxies would be included in the primary 
lens model. 

Figure 4 shows the result for a ,- {{)) as a function of source redshift. It can immediately be seen 
that the contribution of intcrgalactic halos to millilensing can be significant. Previous estimates 
for the amount of substructure required to account for the monochromatic magnification 
anomalies were a few percent. Given that the surface density in the primary lens is k ~ 0.5 
at the critical curve near where the images are located, the siirfacc densities in figure 4 are 
substantial. Just like for Aa{9), aK,{9) will only increase when the clustering of the halos is 
included. 

By comparing figures 4 and 2, it is seen that (Tk(^) has a stronger dependence on smaller 
mass halos, while Aa(9) is heavily weighted toward the high end of the mass range. One can 
also see in figure 4 that the concentration of these halos has an important influence on their 
lensing properties. The redshift dependence is also important. A way to distinguish intcrgalactic 
halos from substructure in the primary lens is to measure how the anomalies depend on source 
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Fig. 6. — The standard deviation in the dimensionless surface density in small scale structure. Everything is the same as in 
figure 4 except that the mass-concentration relation is given by equation (1) without the redshift dependence and with Co = 12. 
The redshift dependence tends to reduce <Tk somewhat. 



redshift. Since the space density of small scale halos becomes larger at higher redshift there 
is a significant increase in a^iz) over the range of observed QSOs lenses ~ 1 to 4). The 

dependence of the normalization of the matter power spectrum may be counter intuitive to 
some. If as is reduced, more of the mass is in smaller halos and since we are only including small 
halos (well below M*, the nonlinear mass scale), is increased. If the current measurements 
of as are taken into consideration and the primordial power spectrum is assumed to be scale 
invariant down to these scales, this effect is relatively small compared with the importance 
of halo concentration. If scale invariance is broken however this could become an important 
dependency. 

A very interesting and illuminating thing occurs when n 7^ 2. Figure 5 shows the same 
quantities as in figure 4 only with the large r slope changed to n = 3. This has the effect of 
cutting the halos off more abruptly and, to conserve mass, increasing their central densities. It 
is clearly seen that ctk is larger in this case. This strong dependence on n demonstrates that it is 
not the inner core of a single halo that is doing most of the lensing, as is implicitly assumed in 
a cross section type calculation. The important property is the profile at radii where the halos 
are likely to overlap on the sky. A cross section type calculation does not account for this since 
there the line of sight is assumed to intersect with at most one halo. 

Figure 6 shows a^lz) with a redshift independent c-M relation. The same c-Af relation was 
used in the simulations of Metcalf (2004). The fiuctuations in k are similar, but larger in this 
case. It is possible that baryonic contraction in the larger halos could convert the inner region 
of their profiles into the SIS form {p{r) oc r~^) at some high redshift. If this is the case, aK,{z) 
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Fig. 7. — The correlation of surface density as a function of angular separation between images. The values for the halo 
parameters and n are listed. The three lowest curves (for most of the 9-range) are for halo mass ranges 10^ M© < M < 10^ Mq 
and the others are for 10^ Mq < M < 10^*^ Mq. The source is taken to be at = 3 and the primary lens is at zi = 0.5. The 
power spectrum normalization is as = 0.9 in all cases. Only one n = 3 curve is included for clarity of the plot. 

will be even larger. If the z = normalization can be established by other means, the redshift 
dependence of the magnification ratio anomalies can be used to measure the redshift dependence 
of the c-M relation. 

As mentioned before, if the fluctuations in k{6) have a coherence length that is much larger 

than the size of the primary lens, they will have no observable effects because of the mass sheet 
degeneracy. Figure 7 shows (JK,{d) as a function of image separation. Most of the fluctuations 
caused by halos of mass < 10^° Mq are on a significantly smaller scale than typical image 
separations. The concentrations of the halos appear to strongly affect only the small scale 
fluctuations below a few tenths of an arcsecond. 

One might be concerned that the variance in k,{0) is not a good quantity for evaluating 
variations. Often, distributions in lensing have long tails which make the variance misleadingly 
large - for example the distribution of microlensing magnification caused by point masses. This 
is not the case here because the halo profile used does not diverge at the center or fall off too 
slowly at large r. This makes cr^ a trustworthy indicator of typical variations. It is possible that 
the halos are highly clustered which would make a^, larger and skew the distribution so that the 
mode lies below the mean. If this were a dramatic effect the standard expression for the angular 
size distance would not apply to most high redshift objects. 

The variance in the shear along random lines of sight, defined as cr^ — (71+72) , can be 
calculated in the same way as tr^ was calculated. However, the requirement that the distribution 
of halo positions be isotropic automatically guarantees that (t^XO) = o''^{0). This is best 
demonstrated in Fourier space so the explanation will be deferred until section 3.2. Despite this. 
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an explicit calculation of (0) has been done with the halo model and it is found to equal a'^(O) 
to very high accuracy which provides a reassuring consistency check. 

The perturbations to the magnification matrix represented in figure 4 are on the several 
percent level. It might seem like this is too small to be of importance, but in fact these 
perturbations can strongly affect the magnifications of the images. To see this, the magnification 
can be expanded around the value for the primary lens alone, m,"^ 

^ ~ 2ni [(1 - ki)6k + 7, • ^7] . (12) 

The variance in the perturbation to the magnification is 



^) ^2M/a«(0)J(l-K02 + ^ (13) 



where the fact that <t'^{9) = (j'^{6) has been used. Galactic sized lenses are well fit by singular 
isothermal ellipsoid (SIE) models, in which case the images form at « 7; « 0.5. In this case 
( ^ ) « ujCTk- For an image with a modest magnification of 5 this will exceed 10% if the 

V / rms 

source is at .Zs ^ 2. For more highly magnified images and/or higher redshift sources this can 
be significantly larger. And if halos have a steeper cutoff than NFW at large r it can be much 
larger. The level of the observed anomalies is on the ^ 10% level (Metcalf & Zhao 2002). And 
direct numerical simulations give similar results (Metcalf 2004). 



3.2. power spectrum approach 

The second order statistics calculated thus far using the halo model can also be calculated 
in terms of the power spectrum of density fluctuations. This may be a more direct, although 
somewhat more obscure, way of interpreting the observations since the halo model contains more 
information than do a-^{0) or Aa{9)'^. The power spectrum also has the advantage of including 
within it the correlations in halo positions which were ignored in the previous section and it can 
be calculated relatively easily from a N-body simulations . With enough data one could think of 
inverting the variance in magnification anomalies to get a constraint on the power spectrum at 
small scales. 

The deflection along a short section of the light path is 

d{0) = 2V_L$ [x{0, z)) 61 (14) 

where ^{x) is the Newtonian potential, V_l is the gradient operator in the two dimensions 

perpendicular to the light path and 51 is the path length. If the imperturbed path is substituted 
in for x{6,z), the potential is Fourier transformed, and all the deflections along the path are 
added up one gets 

a{e) = {SttGpo) [dR{l + z)^^^ f ^6{k) exp (-ik± ■ x±{e) - ikrR) (15) 



^This expansion is not guaranteed to converge when the source is very near to a caustic, but if, roughly speaking, jXiSn < 1 
this should not be a problem. 
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where ^(A;) is the Fourier transform of {p — p)/p, Po is the current density of the universe and R 
is the radial distance. The correlation function for the deflection and any of its derivatives can 
be calculated with the Fourier version of Limber's equation (Kaiser 1992) as is commonly done 
in weak lensing calculations. This is essentially the assumption that the quantities z and D{z) 
change less rapidly with R than 5{x) does. The result is 

Ho \ Us J Jo ^nrn{^ + Zr + Jo « 

(16) 

where P{k, z) is the power spectrum of density fluctuations and Jo[x\ is the zero-th order Bessel 
function. And by taking derivatives of (15), squaring and averaging, one can work out the 
variance of all the parts of the magnification matrix. Specifically 

al{e)=a^^{e) = -4- fdz (l + ^)'^cHt(^,^.)-^ kMky{z,9)]P{k,z,M^^^). (17) 

The first equality comes from the fact that {kf - klfj^ + klkl = {k\ + klfj^ = k'^/A. 

The power spectrum used here should really be a conditional power spectrum. As seen in 
equation (7), the deflection angle is proportional to the mass of the halo so the dispersion can 
be dominated by high mass halos. If no other galaxies are seen near the primary lens or they 
are included in the lens model for the primary lens then there is an upper halo mass (or light) 
cutoff, Mmax7 implicit in this power spectrum. In practice, -Mmax should be a function of z to 
express the detection limits of the observations. 

4. Discussion 

It was demonstrated here that intergalactic matter is expected to cause changes in the 
magnification ratios that are of the same size as those that have been observed. This confirms, 
at least qualitatively, the conclusions derived from numerical simulations (Mctcalf 2004). The 
possibility that the anomalies are from substructures inside the primary lens may not be 
consistent with ACDM (Mao et al. 2004; Amara et al. 2004), although simulations with high 
enough resolution to conclusively rule this out have not yet been done. Probably the best hope 
for distinguishing obscrvationally between these two possibilities is to measure the dependence 
of magnification ratio anomalies on source redshift in a statistical way. 

It was found here that the importance of intergalactic halos to the magnification ratios 
depends on the internal structure of the halos as well as the primordial power spectrum. The 
mass-concentration relation is of particular importance. This relation depends on the formation 
time of the halos (halos that form earlier are denser), the history of mass accretion and, in turn, 
on the slope of the primordial power spectrum and amount of dark energy. The concentrations 
could also change if the dark matter particle has a small mass (warm or tepid dark matter) or if 
it is self-interacting. 

There are several factors that have not been taken into account here and would be expected 
to increase the importance of intergalactic structure. The halos will be strongly clustered which 
will clearly increase (Tk- The scatter in the c-M woiild be expected to increase halos' influence 
on lensing as a population. The contraction of baryons within the halos will increase their 
concentration. Observed lenses (at a higher mass scale) have profiles closer to a SIS profile than 
the NFW profile. SIS halos will produce larger discrepancies in the magnifications. Finally, 
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intergalactic halos will have substructure of their own within them which will increase Gk- For 
these reasons the numbers given in this paper are conservative estimates. 

To go beyond the variance in the image magnifications and to go beyond the unperturbed 
light path approximation numerical simulations are required. It has also been assumed here that 
the source is infinitely small which in many cases will not be a good approximation. One expects 
that when the size of the source is on the order of the size of the halos causing the lensing, the 
variations in the magnification will be smaller. Further simulations of specific lenses are needed. 

Instead of thinking of this effect as lensing by individual galaxies in intergalactic space it 
might be better to think of it as a variation in the angular size distance as a function of position 
on the sky or a variation in the surface density of the universe. A typical line of sight travels 
through many halos (500 with lO'^ Mq < M < 10^ Mq for = 3!). On average they add 
up to the average surface density of the universe and the normal expression for the angular 
size distance calculated with the Robertson- Walker metric is recovered. However, there are 
variations in the surface density at the several percent level and these are enough to change the 
magnification ratios of multiply imaged quasars. This opens up an opportunity to probe the 
small scale, otherwise unobservable, structure in the universe which is shaped by the physics of 
the early universe as well as the history of star formation and reionization. 
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